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Abstract 

Much attention has been paid to search for p-wave superconductivity whose edge 
would support Majorana modes, related to a mysterious type of fermions predicted 
seven decades ago and newly believed to be useful in realizing topological quantum 
computation^^. Possible signs of p-wave superconductivity were observed previ- 
ously in organic superconductors 3 , heavy fermion superconductors 4 , Sr2Ru04 5 , the 
5/2 fractional quantum Hall state of two-dimensional electron gas 6 , and recently in 
hybrid systems composed of s-wave superconductors and strong spin-orbit coupling 
materials^m. However, an agreement has not yet been reached. In this work, we have 
designed and constructed particular types of superconducting quantum interference 
devices at the surface of topological insulator candidate Bi2Te3 by utilizing supercon- 
ducting proximity effect. Two types of Cooper pair interference patterns have been 
observed. One arises from a zero-phase loop inhabited in the bulk of Bi2Te3, and 
the other arises from an arbitrary-phase loop formed involving in the surface states. 
The results are consistent with the scenario that chiral p-wave-like superconductivity 
occurs at the Bi2Te3 surface. 
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Topological insulators (TIs)Pl are a class of insulator with topologically protected con- 
ducting surface due to strong spin-orbit coupling (SOC). The SOC interlocks the momentum 
and spin degrees of freedom of the electrons in the surface states, so that the motion of the 
electrons appears to be helical, acquiring a Berry phase. The helical surface states survive 
in the presence of a conducting bulk, as revealed by both angular resolved photon emission 
spectroscopy studies^ and electron transport measurements^!. 

When the helical electrons pair up via, e.g., superconducting proximity effect, Fu and 
Kane proposed that they will form a chiral p-wave-like superconducting stated Figure 
la illustrates how the pair correlation is delivered from an s-wave superconductor (colored 
in dark blue) to the surface and the bulk of a three-dimensional (3D) TI via the proximity 
effect. While the induced superconductivity in TI keeps more or less to be s- wave-like within 
a certain coherence length from the interface (colored in light blue in Fig. la), at a farther 
distance on the TI surface, however, only those Cooper pairs formed of helical electrons 
survive (colored in orange in Fig. la). The interlock between the momenta and spins of 
the carriers there leads to the formation of a spinless superconducting state whose orbital 
part of the order parameter must be in odd symmetry, presumably chiral p x -Mp y -wave-like. 
Similar to the motion of a single helical electron, a Cooper pair of that state acquires a 
Berry phase while traveling around. 

Experimentally, superconducting proximity effect has been observed between conven- 
tional s-wave superconductors such as Al, Nb, Sn, Pb, NbSe 2 and 3D TIs such as Bi 2 Se 3 
and Bi 2 Te3.^^^1. However, no obvious p x + ip y - wave-like signature has been resolved. To 
identify such a state, one would need to embed a segment of the induced superconductor 
into the loop of a superconducting quantum interference device (SQUID), and then exam- 
ine the interference pattern of the Cooper pairs. This kind of phase-sensitive experiment 
involving a 7r-loop was first proposed for testing the pairing symmetry of heavy fermion 
superconductors^, and was later applied to determine the <i-wave pairing symmetry in 
cuprates^ as well as the possible p-wave pairing symmetry in Sr 2 RuOj^. A 7r-loop can also 
be realized by embedding a quantum dot in the loop and adjusting the even-odd occupancy 
of electrons on the doiPH The most significant feature of a 7r-loop SQUID is its exotic 
magnetic flux dependence of critical supercurrent. Unlike a 0-loop SQUID whose critical su- 
percurrent is maximized in zero magnetic field, the critical supercurrent of a 7r-loop SQUID 
is maximized in finite magnetic fields. 
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FIG. 1: (color online) Schematic illustration of constructing an arbitrary-phase loop by us- 
ing proximity-effect-induced superconductivity on the surface of a topological insulator or related 
material, a, Development of proximity-effect-induced superconductivity from an s-wave supercon- 
ductor (the dark blue disk) to a topological insulator with bulk states. The black arrows represent 
the Cooper pairs. The s-wave character is kept within a certain distance from the interface (the 
regions colored in light blue). Those paired electrons reconciling with the helical surface states (the 
red arrows) at a farther distance form a spinless chiral p x + ip y -wave-like state (colored in orange), 
b, Illustration of a 7r-loop composed of an s-wave half ring and a p x + ip^-wave-like half ring. The 
circular motion of Copper pairs in the latter half acquires a Berry phase of ir. c, Illustration of an 
arbitrary-phase loop, d and e, Two experimentally accessible variations of arbitrary-phase loops 
utilizing the superconducting proximity effect at the surface of a topological insulator, in which 
the two types of segments are cut to small pieces to different extent (multiple segments and many 
small grains, respectively), and mixed up together. 

Figure lb illustrates a 7r-loop composed of a half of s-wave superconductor and a half 
of p x + ipy-wave superconductor. The circulation mode of a Cooper pair in the latter half 
acquires a Berry phase of 7r after turned 180°. To create the latter half, one can make use 
of the strong proximity effect between, for example, Pb and P^Te^ - by using discrete Pb 
segments to bridge the interference path, meanwhile allowing the phase of Cooper pair to 
relax rather than to be fixed to a unitary value of the s-wave state. The idea is illustrated 
in Figs. Id and le. In such rings with mixed s-wave and p x + zp y -wave-like segments, we 
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FIG. 2: (color online) Cooper pair interference in a superconducting quantum interference device 
composed of Pb grains on the surface of a Bi2Te3 single crystalline flake, a, A scanning electron 
microscope image of the device, b, Differential resistance dV/dluas of the device measured at 30 
mK, as functions of bias current /bias and magnetic field B perpendicular to the device plane. 
Two oscillation patterns can be recognized. Pattern A (colored in red) represents the edge of a 
zero-resistance state. Pattern B (with a whitish-blue color) inhabits on a finite resistance state, 
c, The central part of b. The vertical lines are guide to the eyes showing that the peak positions 
of pattern B are horizontally shifted by an amount of 25/2ir ~37% between its positive and the 
negative bias current halves, d and e, dV/dl^^ — -^bias and V — /bias curves, respectively, taken at 
several different fields at 30 mK. 

speculate that the accumulated phase shift of a Cooper pair along a circulation mode is 
in general not necessary tt, but an arbitrary value depending on the turning angle of the 
pair in the p x + ip^-wave-like region. Furthermore, since the clockwise and counterclockwise 
circulation modes are degenerate in energy, spontaneous symmetry breaking happens so that 
the system forms a quantum double-well. 
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Motivated by the above picture, we have fabricated Pb-grain SQUIDs (i.e., the design 
in Fig. le) with different Pb coverage ratios on the surface of Bi 2 Te3 flakes, and studied 
the interference of Cooper pairs via electron transport measurements down to 30 mK in a 
dilution refrigerator. We have also explored the multi-segment design shown in Fig. Id. 
The results can be found in the supplementary material. 

Figure 2a shows the scanning electron microscope image of a Pb-grain SQUID. Between 
the two bulky Pb electrodes there is a ring composed of Pb grains of 10-50 nm in diameter 
and 10-20 nm in height. These grains were formed naturally when Pb was sputtered onto the 
Bi 2 Te3 surface. Figures 2b and 2c show the differential resistance of the device measured at 
30 mK. Two interference patterns can be recognized. One represents a zero-resistance state 
(dV/dluaa < 0.2Q, marked in red), hereafter called pattern A. The other inhabits on a low 
but finite-resistance state, with a whitish-blue color, hereafter referred to as pattern B. The 
periods of the patterns, AB^ = 1.0 ±0.1 G for pattern A (best seen in Fig. 3h) and ABb = 
1.15 ±0.1 G for pattern B, correspond to two slightly different areas of Sa = 21 ± 2/zm 2 and 
Sb = 18.1 ± 1.5/im 2 , respectively, if estimated using O = SAB = h/2e. The oscillations in 
pattern B are relatively regular, with a period which is in good agreement with the effective 
area of the ring^Sl S = 7rDi n D out /4: m 18.8/mi 2 (where Di n = 4/xm and D out = 6/im are the 
inner and outer diameters of the ring respectively). 

The Pb coverage ratio of the device shown in Fig. 2 is 71±2%. By controlling the 
deposition time, devices with different Pb coverage ratios were also obtained and studied. 
It appears that the devices with high Pb coverage ratios (e.g., the ones with continuous Pb 
films- 8 or with a Pb network, see Fig. S2 of the supplementary material) show only pattern 

A, whereas the devices with Pb grains and low coverage ratios most often show only pattern 

B. The results are summarized in Table SI of the supplementary material. 

Pattern A should arise from a 0-loop, since its highest peak is located at zero magnetic 
field in both positively and negatively biased current directions. Compared to pattern A, 
pattern B has the following anomalous features, (i) There is a phase shift between the 
positive and negative current halves, (ii) It inhabits on a finite resistance state, (iii) It 
appears to be more robust than pattern A, i.e, the oscillations are more uniform, with a 
wider envelop against magnetic field and survivable to higher temperatures. Also, as we will 
show later, pattern B survives in larger parallel magnetic fields. Adding to the robustness, 
pattern B was observed on all the six working Pb-grain devices out of seven that were 
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FIG. 3: (color online) a - e, Evolution of the patterns shown in Fig. 2 as a function of parallel 
magnetic fields Bu, measured at 30 mK. The data are plotted according to the color scale in Fig. 
3i. Pattern A is suppressed when Bn > 100 G. Pattern B remains to be visible at Bn ~ 600 G. f, 
The amplitudes of the two patterns as a function of parallel magnetic field at 30 mK. g, A cartoon 
illustrating that a sufficiently large parallel magnetic field will suppress the bulk superconductivity 
(colored in light blue), but leaving the surface superconductivity (colored in orange) less touched, 
h - j, Evolution of the two patterns as a function of temperature at Bu = 0. 

investigated, whereas pattern A was observed only in one of the six. 

Based on the fact that the surface electrons and the bulk electrons manifest distinctively 
in Shubnikov-de Haas oscillation^, we speculate that the induced superconducting states at 
the surface and in the bulk, and the interplay between them, give rise to the two distinctive 
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patterns. The existence of a superconducting surface has been reported by M. Veldhorst et 
alP^. In order to allocate the surface and bulk contributions to the observed two patterns 
in the present experiment, we applied a parallel magnetic field trying to suppress the bulk 
superconductivity in the Bi2Te3 flake. 

Figures 3a to 3e show the evolution of the patterns in several parallel magnetic fields 
B». While pattern A was suppressed when B» > 100 G, pattern B can still be recognized at 
B\\ ~ 600 G. The field dependence of the oscillation amplitudes is summarized in Fig. 3f. 
In Figs. 4b to 4f we further show the data observed on another device, whose pattern B 
survives to ~ 8000 G, a field where the induced bulk superconductivity in TI should have 
been completely suppressed. Figure 3g is a cartoon illustrating that the superconductivity 
in the bulk (colored in light blue) is suppressed in sufficiently large parallel magnetic fields, 
but the superconducting surface (colored in orange) survives. 

Because less diamagnetic energy is needed in field penetration, the effective critical 
field of a thin superconducting plate can be greatly enhanced in parallel magnetic fields^: 
B* = y/2AB c \^/d, where B c and Al are the critical field and the London penetration depth 
of the material in its bulk form, respectively, and d is the thickness of the thin plate. For 
the induced superconducting surface and bulk in Bi2Te3, therefore, the ratio between their 
effective critical fields is B* surf /B* hulk ~ rfbuik/^surf if neglecting any helicity-related differ- 
ence. It roughly agrees with the observations in this experiment, i.e., -Bcsurf being several 
thousand Gauss, B* bulk ~100 G, d hn \ k ~ 100 nm, and c? sur f being the thickness of the surface 
state, which is a few nanometers^. We thus conclude that pattern B is contributed by 
Cooper pairs interference in the superconducting surface. 

In the following we will discuss whether the observed shift of pattern B is an intrinsic 
character of the superconducting surface, or rather caused by the formation of multiple (>3) 
junctions unevenly distributed on the two arms of the ring (e.g., mimic a Josephson junction 
RatchetP^), due to inhomogeneity in the induced superconductivity. 

The fact that pattern B is more robust than pattern A indicates that the induced su- 
perconducting surface is much easier to form an interference loop than the induced super- 
conducting bulk, though they stem from the same Pb grains. It reflects that the Cooper 
pairs on the surface have a longer coherence length than that in the bulk, echoing that the 
unpaired electrons in the helical surface states experience less back-scattering. According to 
our previous studies^ES proximity superconductivity spreads from its parent superconduc- 
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FIG. 4: (color online) Cooper pair interference in a second Pb-grain SQUID, a A scanning electron 
microscope image of the device, b Differential resistance of the device measured at 30 mK. Its 
zero-resistance state (with dV/dl^s < 0.2S7, colored in red) shows only a Fraunhofer diffraction 
pattern but no interference oscillations. Nevertheless, pattern B is still there, demonstrating a 
phase shift of 25/2tt ~69%. c - f , Evolution of pattern B in parallel magnetic fields measured at 
30 mK (only those data within a few Gauss window near B = in b are traced), g - j, Evolution 
of pattern B with varying temperature at B\\ = 0. 



tor to a distance of several micrometers into the TI, a length which is long enough to bridge 
the Pb grains in this experiment to form a uniform superconducting loop. Therefore, for 
those Pb-grain devices whose bulk has already formed a superconducting loop, it is reason- 
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able to believe that their surface states should have also formed a uniform superconducting 
loop, with only two weak links in the middle of each arm. In Fig. S5 of the supplementary 
material, we further show the data obtained on a four-segment SQUID. Despite that the 
four-segment device is almost symmetric as reflected by its unshifted pattern A, an anoma- 
lous pattern B still appears. Therefore, we conclude that pattern B reflects the intrinsic 
property of the superconducting surface. 

At first glance, pattern B in Fig. 4 and Fig. S5 of the supplementary material appears to 
be tilted. However, being shifted or being tilted corresponds to two different mechanisms. 
A tilted interference pattern arises when the applied current exerts a significant amount of 
magnetic flux to the SQUID loopP^. For our devices, however, the flux generated by the 
applied current of the order I c = 1/iA cannot exceed LI c /2 « 0.00100 (where L rs 3 pH is 
the inductance of the ring with a radius of 2.5 /im, and 0o is the flux quanta). Therefore, 
that mechanism of tilting can be ruled out. Actually, all the peaks in pattern B are shifted 
horizontally by a same amount regardless of their height, which can best be seen from the 
data of the four-segment SQUID shown in the supplementary material (Fig. S6). 

A shift in peak position implies that there is an additional phase added to the phase 
quantization condition of the loopP^: 0i + 02 + 27r$/0o + 5 = 2irN, where 0i and 02 are the 
phase differences across the two weak links on the arms, $ is the applied magnetic flux, N is 
an integer, and S is presumably the additional (Berry) phase acquired by the Cooper pairs 
while traveling and turning in the p x + ip y segments of the loop. Since the clockwise and 
counterclockwise circulation modes are energetically degenerate near zero bias current, the 
system is in a double-well state with spontaneously broken time-reversal symmetry mimic in 
an s-wave and <i-wave 7r-loopp^. Applying a bias current tilts the double-well potential via a 
given asymmetric current distribution of the device, thus changes the population/dwelling 
time of the quantum state in the two wells in the presence of inter-well tunneling, leading 
to the seeming tilted pattern (illustrated in Fig. S7 of the supplementary material). 

Another unusual feature of pattern B is that it always inhabits on a finite resistance state. 
Although for each individual device we cannot rule out that the resistance is caused by inho- 
mogeneous superconductivity on or outside the ring, the reproducibility of this feature over 
many devices investigated would indicate that it is a general property of the surface states. 
A finite resistance state at low temperatures is usually caused by phase fluctuations and 
dissipations, especially in superconductors with reduced dimensions. In Josephson junctions 



9 



with a critical supercurrent of 1/iA or smaller, like those in this experiment, phase diffusion 
often happens via macroscopic quantum tunneling, either in a thermally assisted regime or 
in a quantum regime. If phase diffusion is responsible for the finite resistance state of pattern 
B, it should take place between the 2D potential wells in the (0i,02) plane, in a direction 
perpendicular to the clockwise/counterclockwise circulation double wells. In addition to 
phase diffusion, Little Parks mechanism® might also be responsible for the finite-resistance 
state of pattern B if a Berry phase is taken into account, a scenario which is compatible 
with the picture of chiral p-wave-like superconductivity. Further study is still needed to fully 
understand the finite resistance state of pattern B. 

To conclude, we have observed an anomalous Cooper pair interference pattern in addition 
to a normal pattern on our specially designed SQUIDs, and identified that it is contributed 
by the surface superconducting state of Bi 2 Te3. The results are consistent with the sce- 
nario that the proximity-effect-induced superconductivity at the surface is of odd symmetry, 
presumably chiral p-wave-like. The types of the devices in this work might be useful for 
generating and braiding anyonic excitations, and for realizing topological quantum compu- 
tation. 

The method: The Bi2Te3 flakes of ~100 nm thick were exfoliated from high-quality sin- 
gle crystals grown by Bridgman method, and were then transferred onto Si/Si02 substrates. 
Bulk Pb films of ~200 nm thick and grain-like Pb structure were fabricated onto the flakes 
via standard e-beam lithography, magnetron sputtering and lift-off techniques. The electron 
transport measurements were performed at low temperatures down to 30 mK in a dilution 
refrigerator. The differential resistance of the devices was measured as functions of both dc 
bias current and applied magnetic fields in both the perpendicular and parallel directions. 
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1. Material and devices characterizations 

Bi 2 Te 3 single crystals used in this experiment were grown by Bridgman method and were 
proven to be of high quality by X-ray diffraction and electron transport measurements. The 
mobility of the carriers is 2000 - 5000 cm 2 /Vs. Relevent data can be found in Ref. 1. 

After the crystals were mechanically exfoliated into flakes of ~100 nm in thickness, Pb thin 
films (~200 nm thick) and/or Pb grains (10-50 nm in diameter and 10-20 nm in height) were 
deposited to the surface of the flakes via electron beam lithography and sputtering procedures. 
The Pb grains were formed naturally on the Bi 2 Te 3 surface by controlling the deposition time. 
Figure SI shows the optical microscope images of typical Pb-grain devices. 




Figure SI | Optical microscope images of typical Pb-grain SQUIDs. 



2. Pb-grain devices with different Pb coverage ratios 

The central idea of constructing an arbitrary-phase loop using superconducting proximity 
effect is to cut the s-wave superconductor into small segments and mix them up with the Tl 
material to form an interference loop. In order to compromise between the strength of proximity 
effect and the amount of phase shift along the loop, we have tried serveral different Pb coverage 
ratios by varying the Pb sputtering time. The coverage ratio can be determined from the SEM 
pictures of the devices by using softwares such as Photoshop. 

It has to be noted that the sputtering time window for controlling the coverge ratio of Pb in 
our sputtering machine is only ~ 10 seconds, and that there is not a linear relationship between 
the deposition time (hence the nominal thickness of the Pb film) and the coverage ratio yielded - 
it seems that the transition from forming discrete Pb grains to forming a connected netwrok takes 
place rather abruptly. Adding to the difficulties, the morphology of the Pb grains might also 
depend on the surface cleaness of the Bi 2 Te 3 flakes, the vacuum while depositing Pb, etc., 
althought we tried to maintain all the conditions the same throughout the experiment. 
Nevertheless, it would still be helpful to show the statistical data of our investigation. 

Overall we have investigated seven Pb-grain devices in three batches with different Pb 
coverage ratios. Six of them showed pattern B. The situations are summarized in Table SI. 



Table SI | Statistics on devices with different Pb coverage ratio. 
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l c =l- 10 
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l c =0.3- 0.7 uA * 
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l c ~0.5uA* 


no 


no 


Pattern B 


no 


no 


yes 
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* The characteristic currents depend on the detailed geometry of the devices. 



It turns out that once the Pb grains are connected to form a network, which usually happens 
when the sputtering time is longer than ~11 s on a sputtering rate of ~1.7 nm/s in our experiment, 
the interference pattern of the device becomes the same as that of Pb-Bi 2 Te 3 SQUIDs with 
continous Pb film reported in Ref. [1]. No pattern B was observed in this situation even if the 
pattern A of these devices was suppressed in a parallel magnetic field of > 300 G. Figure S2 gives 
an example for such situation. Shortening the sputtering time/lowering the Pb coverage ratio 
reduced the amplitude of pattern A to be below 1 u,A, and allowed pattern B to emerge. With 
further lowering the Pb coverage ratio, only pattern B remained. 
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Figure S2 | Pattern A of a Pb network SQUID (Pb converage ratio >90%) measured at 30 mK and 
in parallel magnetic fields. The zero resistance state is marked in red color (dV/dl<0.2 Q). The 
interference pattern was regular and unshifted/untilted. The pattern was suppressed above ~ 300 
G. No pattern B was observed. 



3. The upper critical fields of Pb thin films and Pb grains 

In the main frame of Fig. S3 we show the resistance of a typical Pb-grain device (the second 
device shown in Fig. 4 of the main manuscript) as a function of magnetic field applied parallel to 
the surface of the Bi 2 Te 3 flake. One can see that the total resistance goes up when the parallel 
magnetic field excesses the upper critical field of the Pb film, ~2000 G, as marked by black arrows. 
And there are another upturns starting at + 1.4 T, indicated by the red arrows, corresponding to 
the lower bound of the upper critical fields of the Pb grains. The inset of Fig. S3 shows the 
resistance transition of a Pb electrode (a continuous Pb film) of that device, whose critical 
temperature is 7.2 K (not shown), and the upper critical field is about 0.2 T at 30 mK, as indicated 
by the black arrows. 

Figure S3 | (main frame) The 
parallel magnetic field 
dependence of the resistance of 
the second device presented in 
the main manuscript, measured 
at 30 mK. Inset, the parallel 
magnetic field dependence of 
the resistance of a continuous Pb 
electrode of that device, 
measured in a 3-probe 
configuration at 30 mK. 




4. Interplay and competition between pattern A and pattern B 



The amplitude of pattern B in this experiment is in sub-uA range. This pattern is observable 
only when pattern A's amplitude is reduced to the same range or even smaller. When the 
amplitudes of the two patterns are comparable, further suppression of the superconducting bulk 
by parallel magnetic field may change the overall ratio of the s-wave segment to the p-wave-like 
segment, thus changing the phase shift in the phase quantization condition of the loop as 
mentioned in the main manuscript. As a result, the relative shift between the positive and 
negative current halves varies within a certain range of parallel magnetic field. Figure S4 shows 
the evolution of phase shift with parallel magnetic fields for the first device presented in the main 
manuscript. 
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Figure S4 | (Upper panel) The relative shift between the positive and negative current halves of 
pattern B as a function of parallel magnetic field marked in the frames, observed at 30 mK on the 
first device in the main manuscript. The lines and the numbers denote the phase shift 28 in unit 
of 2n. (Lower panel) A summary of the phase shift as a function of parallel magnetic field. 



When the amplitue of pattern A is reduced to be comparable with that of the pattern B, its 
shape is influenced by the latter, so that the zero-resistance state takes a zig-zag shape as shown 
in Figs. 3b and 3h of the main manuscript. 



5. Investigation on a four-segment SQUID 

Shown in Fig. S5 are the Cooper pair interference patterns measured from two opposite 
electrodes of a symmetric four-segments SQUID, following the design illustrated in Fig. Id. Similar 
to the Pb-grain SQUIDs, two distinctive patterns were observed. Pattern B, the one with 
whitish-blue color inhabited on a finite resistance state, is oppositely shifted between the upper 
and lower two halves, but pattern A (the red ones representing the zero resistance state) is not. 
The regularly shaped pattern A reflects the uniformity of the superconducting surface involved in 
the interference loop, so that the shifting represents the intrinsic unconventional property of the 
superconducting surface. 
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Figure S5 | Cooper pair 
interference in a symmetric 
four-segments SQUID. Inset, 
Scanning electron microscopic 
image of the device. Main 
frame, differential resistance 
dV/dlbias of the device 
measured at 30 mK, as 
functions of bias current Ibias 
and magnetic field B 
perpendicular to the device 
plane. Similar to the data of 
Pb-grain devices, two 
oscillation patterns can be recognized. Pattern A (colored in red) is asymmetric between the 
positive and negative current directions due to the existence of a hysteretic loop in V-lbias curve. 
Pattern B (in whitish-blue color) appears to be tilted. 
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6. Is pattern B shifted or tilted? 

At a first glance, pattern B in Fig. 4 of the main manuscript and Fig. S5 above appears to be 
tilted. However, this is an illusion. Since the peak height of pattern B in these figures varies 
significatly from the center to the side, it allows us to distingush whether the pattern is tilted, or 
rather being shifted between the positive and negative current halves. Taking the data in Fig. S5 
as an example and replotted them in Fig. S6a, if the interference patterns are tilted in same slope, 
then the peak positions should follow the tilted parallel lines. This is however not true. Obvious 
deviations can be found at the positions pointed by the arrows. In fact, the peak positions keep 
uniformly distributed regardless of the height of each peak, as illustrated with the help of two 
sets (upper and lower) of evenly spaced vertical lines in Fig. S6b.| 




Figure S6 | The interference patterns of a symmetric four-junction SQUID shown in Fig. S5. The 
assistive lines help to show that the pattern B is shifted between its positive and negative current 
directions, rather than being tilted. 



7. Spontanenouse symmetry broken and the formation of a double-well potential for the 
arbitary phase loop constructed on the superconducting surface 



After added a phase shift +8 to the phase quantization condition, the clockwise and 
counterclockwise circulation modes are degenerate in energy, undergo spontaneous symmetry 
breaking to form a double-well system, as illustrated in Fig. S7. The inter-well tunneling allows 
the dwelling time of the system to be varied continuously between the two wells with the applied 
current, via a given asymmetric current distribution of the device, leading to the seeming tilted 
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pattern B. 




Figure S7 | A double-well potential for the clockwise and 
counterclockwise circulation modes of an arbitrary-phase loop. 
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